Incompressible excitonic superfluid of ultracold Bose atoms in an optical lattice: a 
new superfluid phase in the one-component Bose-Hubbard model 



in 
o 
o 

(N 
c3 



o 

B 



Yue Yu 

Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2735, Beijing 100080, China 

(Dated: February 2, 2008) 

We predict that a new superfluid phase, the incompressible excitonic superfluid (IESF), in the 
phase diagram of ultracold Bose atoms in d > 1 dimensional optical lattices, which is caused by the 
spontaneous breaking of the symmetry of translation of the lattice. Within mean field theory, the 
critical temperature of the phase transition from this IESF to the normal fluid (NF) is calculated and 
the triple-critical point of the three phases is determined. We also investigate both configuration 
and gauge field fluctuations and show the IESF state is stable against these fluctuations. We 
expect this IESF phase can be experimentally observed by loading cold Bose atoms into a two- 
dimensional lattice where the atom filling fraction deviates slightly from exact commensurations. 
The signatures distinguishing this IESF from the common atom superfluid (ASF) are that (i) the 
critical temperature of the IEST/NF transition is independent of interaction, unlike the ASF/NF 
transition; (ii) the IESF is incompressible while the ASF is compressible. 

PACS numbers: 03.75.Lm,67.40.-w,39.25.+k 
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Ultracold atoms in optical lattices offer new opportuni- 
ties to study strongly correlated phenomena in a highly 
controllable environment Q, 0, 0, 3 ■ A quantum phase 
transition, the supcrfluid/Mott-insulator transition, was 
demonstrated using 87 Rb atoms 00. Strongly corre- 
lated phenomena for boson systems may be studied theo- 
retically by the Bose-Hubbard model H,E1- Most studies 
focused on the quantum phase transition from the super- 
fluid to the Mott insulator for a commensurate lattice at 
zero temperature. 

Recently, the phase diagram of the ultracold Bose 
atoms has been investigated by Dickerscheid et al by 
a slave boson approach and greatly improved by a slave 
fermion approach p|. It was claimed that the Mott in- 
sulator phase crossovers to the normal fluid (NF) in a 
finite temperature for a commensurate lattice. For an 
incommensurate lattice, it was believed that below a crit- 
ical temperature, the atoms are in the superfluid phase 
even in the strong interaction limit if the normal state is 
thought as a NF. To our knowledge, the phase structure 
of the normal state has not been thoroughly analyzed. 

We define the 'normal' state of the one-component ul- 
tracold bose atoms in d > 1 optical lattices by the van- 
ishing of the order parameter of the common atom su- 
perfluid (ASF). We find that there is an incompressible 
(atom-hole) excitonic superfluid (IESF) phase, due to the 
spontaneous breaking of the symmetry of translation of 
the lattice. For the commensurate filling fraction, this 
IESF phase is right above the zero temperature Mott 
insulator. Moreover, for an integer filling fraction, cre- 
ating an atom-hole pair is accompanied by a double oc- 
cupation. Thus, the IESF phase is enshrouded by the 
crossover from the Mott insulator to NL. To observe this 
IESF phase, the atom filling fraction of the lattice has to 
be incommensurate. As we will see that the critical tem- 
perature of the NF/IESF transition is very low. If the 
filling fraction deviates from the commensuration, say 
one percent from the one atom per site, the critical tem- 



perature of the NF/ASF transition will be higher than 
the critical temperature of the NF/IESF transition for a 
quite strong interacting strength. Thus, to observe such 
an IESF state, the filling fraction has to be not exact but 
very close to an integer. The justification to distinguish 
this IESF from the ASF is that the critical temperature 
of the IESF/NF transition is not interacting-dependent 
while the critical temperature of the ASF/NF transition 
decreases as the interaction strength increases. On the 
other hand, the IESF is incompressible while the ASF is 
compressible. 

We investigate the ultracold atoms by the Bose- 
Hubbard model 0, |(| with Hamiltonian H = 

i°j _ ^E"i+ lEi^K - 1) +Vtra P - Here a] 
is a Bose atom creating operator on site i and the sym- 
bol (ij) denotes the sum over all nearest neighbor sites. 
/i is the chemical potential, t is the hopping amplitude 
and U the on-site interaction, which are determined by 
the optical lattice parameters and the s-wave scattering 
length of the atoms. We study d > 1 optical lattices only 
in this work. We first consider a homogeneous case with 
V t 



rap 



0. In order to explore the finite temperature 
behavior, we decompose the boson operator by the slave 
fermions, a\ — J^^Lo v / aH-Tc^ +1 jC a) j, where the slave 

fermion operators c a ,i obey {c Q .i,cjg } = 5 a p6ij . As 
the auxiliary particles, they have to obey the constraint 

J2 a n i = E Q c a i c <*,i = 1 on each site. In the slave 
fermion language, the partition function of the_system 
reads Z = Tte~ pH = / Dc a Dc^ a DX e" SB where " 

f 1/T f 

S E [c a ,c a ,X}= dr < c a ■ [d T — fx 



U 



a(a - 1) + iXi]c at i - i ^ A, 



(1) 



(ij) a/3 



P,j C P+l,3 
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where Xi is a Lagrange multiplier field and n is the atom 
filling fraction of the lattice. 

To study the ASF/NL transition, we decouple the 
four slave fermion term by introducing a Hubbard- 
Stratonovich field which is a bosonic field and may 
be identified as the order parameter of the common atom 
superfluid. In our recent work, we have given the details 
to describe this phase transition by the slave fermion ap- 
proach |j| • We find that using a finite type slave fermion 
approximation, say the maximal qlm = 6, the well-known 
special points in the phase diagram may be well repro- 
duced |9|- For example, for the filling fraction n = 1, the 
non-interacting Bose gas in a three-dimensional lattice 
has the critical temperature T c i=s 7.08t while our calcu- 
lation for olm = 6 gives T c « 7.01i; the zero temperature 
critical interacting strength from the superfluid to Mott 
insulator is U c /zt 5.82 in the mean field theory [Hj 
while our result is U c /zt « 5.9. 

We now focus on the 'normal' state in which = 0. 
We decompose the four slave fermion term by auxiliary 
Hubbard- Stratonovich fields Xap,ij and Va/3,ij, namely, 
the Lagrangian of the system is rewritten as 

L = y^y^<4jd r - «M + %a{a - 1) - i\]c m 

i a 
i (ij) a/3 

- f)ip,ij(Xa+w+i,ij ~ c I+i,i c /m,i) + h.c], (2) 



where t a @ — t^Ja + IV f3 + 1. Integrating over 77 and 
X, the above Lagrangian restores the original one. The 
auxiliary fields Xa/3,ij, Vap,ij an d Xi may be rewritten 
as X a 0,ij = Xaf3,ije iAii , T] a 0,ij = Vap,ije lAi:i and A, = 
X+Aio- Consider the mean field state, Xaf3,ij ~ Xa,ijSa/3i 
Va(3,ij ~ Va,ijSa(3, and Xi ~ X. The mean field Lagrangian 
reads 

L M f = C L M « C ™ - H H \\t<*(rfx,H ~ Xo,y)Xa+i,i 



(3) 



where M Q = 9 r — afi+^-a(a— 1)— iA. The free single slave 
fermion Green's functions can be read out from @,i.e., 
D a (x,T) =TJ2 n e IU " T -D a K) with the matrix 



D a (uj n ) = [(iu n + —a(a - 1) - fj, + iX)5ij 

+ t a -ivl-i, tJ +t a V a +i,ji}~ 1 - (4) 



The mean field equations are given by Xa,ij ~ Vaij and 

Va,tj = (cL c aj) = Tj^n^ij^i^n)- Near the critical 
temperature, these mean field equations have the follow- 
ing solutions 



?7Q ' y C ~ ( e /3c a) (f a(a-l)-afj,-iX) + j\2 ' 



(5) 



For a > 2, the critical temperatures vanish, T, 
while 



.(«) 



77 f 

T (o) _ Jh_ f 

C ~ A ' 

4% 



T(l) « 



(6) 



if Pc/j, 3> 1. The coupling between 770 and 771 means 
771 = 2jyJ and T c (0) = T c (1) = T c = t/2 HJ. Changing 
the variable Ao = —ir/o/t and using 771 = 2?7q, the mean 
field free energy near the critical temperature may be 
expanded by I2| 



F = ^(2t- 1 A^.-r- 1 |A , J | 2 ) 



(ij) 



7f^ E A 0i ij AojfeAo^jAo,, 
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(7) 



The first term in (JZJ) implies only a real Ao,y minimizes 
the free energy, which is given by 



F = -Mt 2 SZ/S 4 



(8) 



where r 



T c -T . 



£2 and 5*4 are the numbers of the non- 



zero terms of the first and second summations in Q ■ Due 
to the filling fraction very close to n = 1, we consider all 
holes in the lattice are isolated. There are two solutions. 
One is the equal-bond state. Namely, for a hole at site 
i, A 0)iii±Ta . = -Ao t i : i± Ty = A for all nearest neighbor 
sites j and all others are zero. Another solution is an 
atom-hole exciton state and Aq^- = A is not zero only 
if a link (ij) is occupied by an exciton, which breaks the 
symmetry of translation of the lattice in the link direction. 
The free energies per site for the equal-bond and exciton 
states are degenerate, e.g., for a two-dimensional square 
lattice with a low hole density, 



feb = fe = -6n°tT 2 



(9) 



where n° is the hole density. There is, however, a config- 
uration fluctuation to the equal-bond state. When two 
holes are close so that they are in the diagonal line of 
a plaquette, the equal-bond state will gain the free en- 
ergy an amount |tr 2 while the free energy of the exciton 
state does not change. Multi-hole configurations further 
raise the energy of the equal-bond state. On the other 
hand, due to the gauge fluctuation as we shall see the 
equal-bond state may not be stable. 

To study the gauge fluctuations, we first calculate 
the dispersion relations of the slave fermions for these 
two mean field states. For a two-dimensional square 
lattice, we divide the lattice into two sublattices, even 



and odd. The 
noted by either 
even or odd 
poncnts e Q; fe 
1 



slave fermion operators c a ^ is 



de- 



respectively. 
1 T- 



;, corresponding to i 
Using the Fourier com- 
e~ tk ' R 'e a i and d a u = 



E 



i£odd 



the hopping term of the Hamil- 



tonian in the mean field theory may be diagonalized 



H t = 



E 

a=0,l;fc 



|A (Q) (fc)|(C a , fe C a , fe -£ Q , fe W), (10) 
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where A<°)(fc) = Ai(fc) and A«(fc) = A (fc). A a (k) = 
A aA e lk * - A a , 2 e- lk y + A ai3 e" ife * - A aA e ik v where 1,2,3 
and 4 denote four adjacent sites around a hole. The 
diagonalizcd operators are defined by 



»A 



ot.k 



2 ( ^ fc+ |A a , 



Ca,/s); 



«A 



|A a , 



(ii) 



For the equal-bond state, the slave fermion dispersions 
are e e fc, Q (fc) = ±2|A a |(cosfc a; +cos£; y ), which have gapless 
excitations. For the exciton, the slave quasi-fcrmions are 
dispersionless with ed,a{k) = ±|A Q |, which means that 
exciting a slave quasi- fermion has an energy gap |A a |. 
However, the slave fcrmions are auxiliary particles and 
the single slave fermion Green's function is not gauge 
invariant. To see the real quasiparticle excitation, one 
has to calculate a gauge invariant Green's function. The 
single-atom thermodynamic Green's function is gauge in- 
variant, which reads 



R v (a + l)(n*-n^) 

K ' Z -w-a[/+|A(«+ 1 )|-|A(«)|+ M + zO+ 



(13) 



where the order parameters A^ ) = Ai, A^ = Ao and 
vanish for others. 

Thus, the lowest energy quasiparticle excitation at fi- 
nite temperature is 



s(k) = |Ai(fc)| - |A„(fc)| = |A (fc)|. 



(14) 



This implies the low lying excitations of the exciton state 
has a gap |A |, which is equal to t at T = , while 
it is gapless for the equal-bond state. The second level 
excitation spends an energy U — |Ai(fc)|. 

We now discuss the gauge fluctuations. The gauge fluc- 
tuations come from the gauge field Aij and Aoi- For the 
equal-bond state, the gauge fluctuations may be more 
serious than in the t-J model because we do not have 
a large N limit. However, the existence of the gap in 



N 




FIG. 1: The phase diagram for 2d for n=l. ASF, NF, IESF 
and MI are standing for the atom superfluid, normal fluid, 
incompressible excitonic superfluid and Mott insulator, re- 
spectively. The shadowed area is the crossover regime from 
MI to NF. 



(T(a(x, 4)0+ (0, 0))) = V^TlV/3 + 1 

a/3 

(T(ci(x,t)c a+1 (x,t)cl +1 (0,0)cfs(0,0))). (12) 

If we consider the lowest lying excitation, only relevant 
propagating processes are those which do not raise an 
additional on-site energy U. That is, at t = 0, positions x 
and have the occupation numbers a and a + 1, respec- 
tively, and after a time t, one atom propagates from to 
x, the occupation numbers become a + 1 at x and a at 
0. The approximation we are using yields the following 
factorization of the four slave fermion Green's function, 
£ a/3 V^TlV?TT(T( C t ( x , t)c a+1 (x, r)4 +1 (0, 0) C/3 (0, 0))) 

£> + l)<T( C t(x,t)c a (0,0)))<T(c a+1 ( a;) t)4 +1 (0,0))). 
The corresponding retarded Green's function is given by 
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FIG. 2: the phase diagram for 2d for n=0.999. The inset is 
the value of U c in the triple-critical varying as the atom filling 
fraction. 



the exciton state strongly suppresses the gauge fluctua- 
tions. Integral over the slave fermion field, one can get 
the effective action of the gauge field 

S[a] = — J d 2 ka^(k,iv)a u (k,u;)n^(k,uj), 

n a (p + k/2)-n a (p-k/2) 



+ 



d 2 p 
2 e 



a 

^ de a de a 

dp„ ' dp^dp v 



w-[e a (p + fc/2)-e a (p-fc/2)] 



-n(p) 



(15) 



where a M (fc,w) is the Fourier component of the con- 
tinuum limit of the gauge field. The dispersionless of 
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the exciton state means II M „ = which suppresses 
the gauge fluctuation. The physical meaning of 11^ is 
the atom density-density and current-current correlation 
functions. The vanishing of the density-density correla- 
tion function yields the incompressibility of the state. 
For the equal-bond state, the gauge field propaga- 

tor n^ 1 = ( a /i(«z) a f (-?)) = (<W - %q. v /q 2 )D T {q) with 

D T (q) ~ [w/q- Xdq 2 ]" 1 for u> < q The gauge 

fluctuation may renormalize the hopping amplitude t to 

t = t(e tAiJ ) = te~( Ai ^ and the critical temperature is 
reduced to T c = t/2. 

To sum up, both configuration and gauge fluctuations 
destroy the equal-bond state while the exciton state is 
safe. The excitons form an incompressible fluid. Fur- 
thermore, because Ai = 2Aq ^ 0, the order parameter 
{a\aj) of the exciton condensate is not zero. This means 
the excitons are condensed when T < T c . This is an 
incompressible excitonic superfluid (IESF). 

The phase diagrams of the system now may be depicted 
according to the above discussions. Figgis the phase di- 
agram for the filling fraction n = 1 for two-dimensions. 
There is a triple-critical point with U c /zt w 5.8 and 
T c /zt = l/2z. The IESF is right above the Mott in- 
sulator. Therefore, this IESF may not be observed in 
the commensurate filling fraction because the IESF is 
still in the range of the practical Mott insulator. The 
crossover from the Mott insulator to the normal fluid en- 
shrouds this IESF phase. To observe the IESF phase, one 
should work in an incommensurate filling fraction. How- 
ever, the critical temperature T c = t/2 is so low that the 
U c / zt <~ 15 for n — 0.99 in two-dimensions. In the inset 
of Fig. we show the critical U c for 0.99 < n < 0.999 
. Thus, the IESF can only be observed when the filling 
fraction deviates slightly from the commensuration. In 
Fig. [3J we show the phase diagram for n = 0.999 for 
two dimensions. For three dimensions, the triple-critical 



point is U c /zt ~ 16 for n = 0.999, much stronger than 
that in two dimensions. 

We now briefly discuss the inhomogeneous case with a 
trap potential Vt rap — VYlj r i n i where is the distance 
of an atom from the center of the trap. In the regime 
of U we are concerning, there is a central Mott plateau 
which is surrounded by a superfluid ring Our re- 

sult points out that there, in fact, are two superfluid 
rings. The IESF ring right surrounds the Mott plateau 
and the ASF ring is outside of the IESF ring. As U is en- 
hanced, while the radius of the Mott plateau almost will 
not vary, the radius of the zero compressibility valley will 
be enlarged because the IESF ring becomes wider. This 
requires an examination by a quantum Monte Carlo cal- 
culation, which will be reported elsewhere because it will 
take a long computing time period and more space to 
describe. 

In conclusion, we predicted a new superfluid phase, 
the incompressible excitonic superfluid, in the one- 
component Bose-Hubbard model for d > 1 if the on-site 
interaction is large enough. This may be observed when 
the atom filling fraction deviates slightly from a commen- 
surate one. Two characteristics of the IESF may be dis- 
tinguished from the atom superfluid: (i) the incompress- 
ibility of the IESF with a gap - O.l^K at T = for 87 Rb 
(The Mott gap > 2.5/iK.); and (ii) the [/-independence 
of the critical temperature T c — t/2 ~ 0.05/xK. We calcu- 
lated the free energy in a mean field theory and analyzed 
the configuration and quantum fluctuations. We intend 
to examine this IESF phase by a quantum Monte Carlo 
calculation. 
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